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Abstract. We show that non-relativistic Quantum Mechanics can be faithfully represented 
in terms of a classical diffusion process endowed with a gauge symmetry of group Z4. The 
representation is based on a quantization condition for the realized action along paths. A lat- 
tice regularization is introduced to make rigorous sense of the construction and then removed. 
Quantum mechanics is recovered in the continuum limit and the full U{1) gauge group sym- 
metry of electro-magnetism appears. Anti-particle representations emerge naturally, albeit 
the context is non-relativistic. Quantum density matrices are obtained by averaging classical 
probability distributions over phase-action variables. We find that quantum conditioning can 
be described in classical terms but not through the standard notion of sub cr— algebras. Del- 
icate restrictions arise by the constraint that we are only interested in the algebra of gauge 
invariant random variables. We conclude that Quantum Mechanics is equivalent to a theory 
of gauge invariant classical stochastic processes we call Stochastic Mechanics. 



The question of whether the Schrodinger equation can be interpreted as a classical diffusion 
attracted much attention since the discovery of quantum mechanics. Semiclassical expansions 
are introduced in (Weyl 1927) and (Wigner 1932) while a Hilbert space approach to quantum 
mechanics is introduced in (Koopman 1931) and (von Neumann 1932). These approaches are 
useful to bridge the gap between the classical and quantum formalisms but do not establish 
an equivalence. A more radical departure is in (Bohm 1952) and (Fenyes 1952), where new 
models are introduced which are not entirely consistent with standard Quantum Mechanics. 
This line of research was then greatly expanded upon in (Nelson 1967) in what became known 
as Nelson's Stochastic Mechanics. See also (de la Pena-Auerbach 1970), (Jammer 1974), (Guerra 
and Ruggiero 1973), (Bacciagaluppi 2005). The key difficulty stems from the inherent differences 
between Quantum and Classical Probability also revealed by Bell's inequalities regarding hidden 
variable theories, see (Bell 1966). 

In this paper we introduce a purely classical representation for quantum mechanics which is 
entirely faithful. We do so by considering non-relativistic quantum mechanics but suspect that 
the construction is of general validity. The classical diffusion on which we base the analysis 
has a Z4 gauge symmetry and it is this symmetry which is responsible of the subtle differences 
between Quantum and Classical Probability theory. On the quantum side, this gauge symmetry 
is also related to the U{1) gauge symmetry of electromagnetism. In Section 1 we summarize our 
results and in the following sections we give details. 
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1. Quantization Condition 

Consider a spinless particle of mass m and charge e in an external electromagnetic potential 
{ A (x) , (j){x)) . To quantize the motion, consider the following diffusion equation: 

e — > /fiT 

(1.1) dlct^ -Arx)dt+\—dW. 

cm V m 

The action realized along a path is formally given by the process 

(1.2) ^* " ? (^) '^^ c^^^^ ■ " 

The trouble with this equation is that, since the paths of the Wiener process are rough, 
the realized action is infinite. The singularities originate from the kinetic term only, while the 
other two terms are mathematically well defined as stochastic integrals. A regularization is thus 
required. 



By using equation (1.1) we can rearrange this expression as follows to extract the first singular 
term: 

where 

e2 _> 

(1.4) V{x)^ -^A{xfdt + e^{-^t)- 

To render all expressions finite and extract the singularities, we fix an elementary space scale 
a > and an elementary time interval 5t > Q. The strategy we follow is to discretize the diffusion 



process in (1.1 1 so that it evolves on {a'L)'^. Then, we regularize also the definition of kinetic 



terms in the action by setting 

(1.5) s,-t f( -"-'--^" + - [ 

As a criterion for the regularization scheme, we ask that the action be quantized according to 

(1.6) St^So^ fvixt)dt + hnt + o(^-^) 

Jo \a m J 

where nt is an integer valued process. In Section 2, we show that it is possible to achieve this 
objective as long as one chooses the time discretization interval to be 

(1.7) H^"^. 

Next, we notice that the joint process {~xt,nt) is translation invariant in the nt direction. 
This derives from the fact that the action has the form of an integral extended over a path: as 
time advances, the realized action is updated on the basis of the most recent changes and without 
memory effects. In (Albanese 2007&) and (Albanese 2007a) processes with similar symmetries 
are called Abelian as they are associated to a commutative operator algebra. Finally, we notice 
that the dynamics of rit can further be restricted to Z4 by using mod 4 periodicity without 
compromising the Abelian character of the dynamics. This gives rise to a gauge symmetry for 
the classical process {x t, nt). Details are in Section 3, but we anticipate here some conclusions. 

The joint kernel uCa?, 0; 1?', n; i) obviously preserves probability in the classical sense. How- 
ever, as consequence of the Abelian symmetry of the action integral. Quantum Mechanics makes 
a wondrous appearance through the following equation: 



(1.8) e*"(^, Ic') = exp ( (1 - i) ^ + 2Kot ) ^ ^"u(^, 0; n; t). 

n=0 
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On the right of this equation we see the classical probability distribution function for the joint 
process. Here, Kq is an energy threshold we define more precisely below. On the left hand side, 
H is the quantum mechanical propagator for the Schrodinger operator 

Z?'^ — > — > 

(1.9) m=-—Aa+—A{^)-Wa + V{x). 

Zm cm 

This formula depends on the lattice spacing a and the argument of the exponential factor in 



equation (1.8 1 diverges as a — > 0. However, the resulting Hamiltonian is well defined, its limit is 
regular and converges to the usual Hamiltonian for non-relativistic Quantum Mechanics. 

These equations are intriguing as they unveil a mathematical relationship, but one more step 
is needed to arrive to a proper and physically grounded representation of Quantum Mechanics. 
Again, the existence of the Z4 gauge symmetry motivates us to revise the starting point and 
complicate a bit the classical model by creating two copies of it. As we explain in Section 4, we 
introduce a second independent joint process (^tjt't). Here satisfies also an equation of the 



form ( 1 . 1 1 except that it is driven by an independent Wiener process and is defined similarly 
to fit except that it corresponds to the time-reversed process. We then look at the classical 
stochastic process with probability distribution function Pc{S,t, i^t,xt,nt',t). Our main result can 
be expressed as follows: 

Theorem 1. The operator of matrix 

(1.10) p,(6, xut)=exp(2 -+ AKot V i""-" pd^ut^t, a;*, n*; t) 

\ ma'' I ^-^ 

^ ^ n.u 

solves the Quantum Mechanical equation for density matrices 

(1.11) pq{t) = e-™p(0)e™. 



Furthermore, all quantum mechanical density matrices can he written in the form (1.10) by 
averaging some classical joint probability density pc{^t,i^tTXt,nt',t). 

The proposed representation is thus completely equivalent to ordinary non-relativistic Quan- 
tum Mechanics. We apparently accomplished the feat of obtaining Quantum Mechanics out of a 
limiting procedure involving two independent copies of a classical stochastic process for two pairs 
of joint diffusions of a particle and the corresponding realized action on the particle world path. 
However, we cannot jump to the conclusion that there is a one-to-one correspondence between 
classical and quantum events. There are subtleties as the classical dynamics of quantized action 
is subject to a Z4 gauge symmetry. This symmetry is of pivotal importance in the construction 
and develops into the full fiedged U{1) gauge symmetry of quantum electro-dynamics in the 
quantum representation. It also affects the base mathematical structures and definition of the 
algebras of quantum events. A measurement apparatus built with physical matter is inescapably 
subject to the same Z4 gauge symmetry that we leveraged on. Hence, physical observables which 
are measurable ought to be gauge invariant. One thus needs to consider the cr-algebra of classi- 
cal events along with the algebra of random variables with are gauge invariant. This is a subtle 
restriction that invalidates the standard constructs in Classical Probability such as that of con- 
ditional probability. This is also at the origin of quantum coherence, quantum entanglement 
and other departures of Quantum Mechanics from the standard Classical Probability. But when 
gauge symmetries are accounted for, faithful mathematical equivalence results. 

The classical representation is not only mathematically complete and faithful, it also contains 
a few extra bits of known Physics. By changing the factor i" into i^" in the right hand side of 



equation (1.8 1, we find another interesting equation 



(1.12) e'"'(^,^') =exp ( {1 - i)-^ +2Kot]y" i-"uCx,0;lc' ,n;t). 

\ ma-' I ^-^ 



n=0 



where H' is the PT reversal of H, i.e. the operator obtained by inverting both time and space 
coordinates. Equivalently, it is the charge conjugate operator associated to antiparticles, as seen 
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in relativistic quantum theory. Similarly one can conceive of seas filled of particles and anti- 
particles along the same lines. We conclude that relativity is not responsible for the existence of 
antimatter, gauge symmetries are. 

2. Lattice Regularization 

Consider physical space discretized on the lattice (aZ)^ and consider the process described 

by the Markov generator 

(2.1) £(^5 = 2^A„(^; + ^^-(^) • V^(^; + ^^^(^) • V-(^; 
where 

/r. r.x A -^/\ K'^' - a? - aTci) + 5{'x' - Ic + aet) - 2S{'x' - Ic) 

(2.2) Aa{x;x) = }^ -2 



(2.3) 

We set 
(2.4) 



Vl+(^;^') = 

vr(i?;i?') = 



S^lc' — Ic — a~e i) — S{~x' — 'x') 
a 

5{lc' — 'x+ aei) — 5{^' — Ic) 



Af{'x)=ina,x{A%x),i)) and A- {~x) = inax{- A\x),i)) 
Let us rearrange this formula as follows: 



jC{'x; 'x') = 



2maP' 



E 

i=l 



(2.5) 



+ ^1 + ^A+(^)^ 5{'x' - "x + aei) 
_(2 + :^ArP) + ^^A-(l?))*(a"-l?) 



Here, 5 is the function such that ^(0) = 1 and 5{^) = if a; ^ 0. 

Next, let us introduce an integer valued process rit such that the joint process {xt,nt) is 
defined by the following lifted generator: 
^ 3 



C{'x,n; lc',n') = 



E 



2(ie 

(5(n' - n + 1) + — A~ (■?)5(n' - n) ) bi^' -Ic - aei) 



+ (^5{n' -n+l)+ '^Al{lc)S{n' - n)^S{lc' -Ic + a~ei) 

-(2+ ^A-(^) + ^4^^)) Si^' - ■^)S{n' - n) 
+ ^V{'x)5{'x' - 'x){5{n' - n - 1) - 5{n' - n + 1)) 
+ ^Ko{5{n' - n - 1) + 5{n' - n + 1) - 25{n' - n))5{'x' - Ic). 

(2.6) 

Here. 5 is the function such that (5(0) = 1 and 5{^) = if a; ^ 0. 5{n) = 1 if n = and zero 
otherwise. Furthermore, we assume that 

(2.7) |y(^)|<2i^o and i^o < 4^- 

The first condition ensures that off-diagonal elements in the Markov generator stay positive and 
the upper bound on Kq is required in an estimate below. As a J, 0, the energy cutoff Kq also 
diverges. 
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The elementary propagator over a time interval St is given by 
(2.8) "stCa?, n; 'x', n') = e^*'--{~x, n; ~x', n'). 

We find 

Et [nt+5t - nt\nt = n,!? i = Ir] = ^ usti^E ,n; ~x\n'){n' - n) 

(2.9) 



3hSt V('x)St ^fst^h'^ 

^ h O' 



where we made use of the bound in (2.7 1. We also have that 



t+st 



St+st 



(2.10) 

Hence, if 



5t 


m 1 








.2" I 



5t 



+ —A{-^) \ -y(x) 

cm ' 



St 



= — ^ ust{~x .n^lc' ,n') 



3 v{~x)5t ^fst^n^ 



nt = n, xt = X 



ml X — X 

T V Jt 



+ — A(^) -y(l?) 



cm 



St 



(2.11) 



St 



then the regularized version of the action is quantized in multiples of h in the sense of equation 

I 2 

(|1.6|). If the limit as a | is taken while holding the ratio |^ fixed, the action will stay correctly 



quantized. 



3. The Joint Process for the Position and the Realized Action 



Notice that the lifted generator in (2.6 1 is defined in such a way that the Markov generator 
is invariant under translations in the direction of n. Using the terminology in (Albanese 2007&) 
and (Albanese 2007a), the pair (l?i, St) is an example of Abelian process. We can thus single 
out a sector with respect to the translation symmetry in the n direction. 

Before proceeding, let us also notice that the lifted generator can be interpreted as describing 
a dynamics on the reduced configuration space (aX)'^ x Z4 by identifying values of n which are 
equal modulo 4. As far as this generator is concerned, we could even restrict n to Z3 but then 
we would not recover Quantum Mechanics. The additional symmetry that Z4 has with respect 
to Z3 appears to be essential in the argument below. 

Consider the partial Fourier transform operator of kernel 



(3.1) 



.F(^,p;^',n) = e-''f"4 



where p G [0,27r). Partial Fourier transforms in the n variable are a block-diagonalizing trans- 
formation for the lifted generator. Let us introduce the operator £ such that 



(3.2) 



{TLT ^)(^,p; x",p') = C{'x,^';p)Spp>, 
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I.e. 



2ma? 



E 



2ae 



-ip 



2)5{~x' -Ic). 



(3.3) 



The sector with p = ^ is special because in this case we recover the quantum mechanics Hamil- 
tonian, i.e. 



C\ X ^ X , ^ 



n 



(1-z) 



dh 



2Ko 



S{^' -Ic). 



fi^ ie/i — ^ — > 

— Aa(^;^') + — M~^) ■ a{~x\lc') + V{x)5{~x' ~lc). 
Am cm 



(3.4) 
where 

(3.5) H(^;^') 

Ain cin 

Let US notice that the kernel of the stochastic process on the principal bundle (aZ)^ x Z4 is 
given by 

(3.6) u(i) = e*^ = exp {tTCT-^)T = T-^e^^T. 
Similarly, we have that 

(3.7) 

dh^t 



7/1 X ^ X ' 2 ' ^ 



exp \tL\— \ I ("a?, ~a?') = "S^ i^uCx, 0; ~x \ n; t) — exp f iffl — (1 — i) ^ ) 



More explicitly, the quantum mechanical kernel can be reconstructed from the probabilistic 
kernel as follows: 

00 



(3.8) 



im{ 



dh t 

{~x,^') = exp ( (l-i) ^+2KQt] Vi"e*^(x',0; x",n;t) 



ma 



This is an intriguing formula as it relates a quantum mechanical observable to a classical diffusion 
kernel. 



4. Density Matrices 

We have seen that because of the Abelian character of the process ('x'f , ), i.e. of the trans- 
lation invariance with respect to the n coordinate, the Fourier momentum p conjugate to n 
is conserved. It is natural to make the hypothesis that only observables whose expectation is 
indifferent to the action of the Z4 gauge group are physically measurable. In fact, a measure- 
ment apparatus itself would have to be a physical system subject to the same symmetry. The 
problem with this loose statement is that one cannot make it mathematically precise if one uses 
a representation with only one classical particle. To make things work, we need to complicate 
the construction and assign to each single particle in the quantum representation two particles 
which evolves independently in the classical representation. 

As a full classical description of a quantum particle we take a quadruplet ( ^ t, ^"4; ~a?t, nj). 
Each pair {~xt,nt) and {£,t,i^t) is postulated to evolve independently. The first pair evolves 
according to the generator in (2.6). For the second pair, we introduce a conjugate process whose 
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dynamics also provides a different lifting of the same base process in ( 1.1 ), given by the following 
generator: 

3 



E 

i=l 



Onp > 

6(1.' - - 1) + —Aii e )S{^' _ ^) )5( e ' - e - alt,) 



2ae 



+ ^nmt' ~imi^' - ^ + 1) - s{^' - ^ - 1)) 



(4.1) 

The Fourier transformed kernel is 



+ (e^^+^A+(0)^(e'- e +a^.) 



:ifo(e 



(4.2) 



Notice that this Fourier transformed generator is the complex conjugate of the generator in ( 3.3 1 . 
Since the Hamiltonian is self-adjoint, we have that 



(4.3) 



where EI is the Hamiltonian operator in (3.51. 

Consider the joint classical probability density for the process ( ^ t, I't', ~x ti nt) and write it as 
follows: 

(4.4) pc{^ ,i^;^,n;t). 

Also form the following reduced matrix 



(4.5) 



p,(t,^;i)=exp(2 ^+AKot) V ^"-^(T,'^ 

V ma'' I ^-^ 



~x , n; t). 



By the results in the previous section, we know that 



(4.6) 



-ita 



p{0)e' 



itm 



This means that the operator p(t) is a properly defined quantum mechanical density matrix. 
We thus proved that a quantum mechanical density matrix is given by a phase average of the 
classical diffusion times an exponential factor of time. 

The density matrix we constructed may or may not correspond to a pure state. The derivation 
holds in general. A pure state is obtained if there exists a wave-function ijj{x, t) such that 



(4.7) 



Pq{ C ,^;i) = '0(C,i)*V'(a:,t). 
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A generic classical joint density will give a valid generic quantum matrix once phase averaging 
is carried out. 

A physical observable in the classical description is given by a random variable of the form 

(4.8) F.Ct, v; -i, n) = i^-^F.il^- 1). 
Let us denote this correspondence with 

(4.9) 7:F, = F,. 

In this case, the classical expectation of the classical observable coincides with the quantum 
expectation, i.e. 

(4.10) ^pg(e';^)-F9(T;^) =exp(2^+4ifoi') XI Pdl ,n){TTFq)(^ ,n). 

Classical dynamics is recovered in the limit as the mass m ^ oo or as /i ^ 0. In this 
limit, phase averaging speed becomes infinitely fast and the quantum density matrix pg{ ^ ; "x ) 
is localized along the line x = y. In this case, the standard stationary phase argument shows 
that trajectories minimize the classical action. 

5. Classical and Quantum Conditioning 

In classical probability one follows Kolmogorov in giving a structure of cr-algebra to event 

space. The situation here is rather simple as we arc making use of a space discretization and 
could safely assume also a finite volume cutoff. So one can safely say that the set of all possible 
events is given by the set of all finite sets of quadruples (^, v, x, n) G {oLY x Z4 x (aZ)** x Z4 
contained within a given large box. The point is that the set of all random variables on is too 
large as not all functions correspond to physical observables. What we need is to restrict random 
variables to only those of the form Fq = -kF^ where Fc{^,x) is a function independent of v and 
n. The problematic part of this is that this reduction cannot be accomplished by reducing the 
sigma-algebra to a quantum subalgebra and declaring that quantum observables are the ones 
measurable with respect to the sub-algebra. This standard way of phrasing conditioning fails 
here as there is no sub-algebra such that all the corresponding measurable ftmctions are precisely 
those of the form Fq = ttFc while non-gauge invariant functions are non-measurable. This gauge 
symmetry is the reason why quantum conditioning is not equivalent to classical conditioning. 

When conditioning to a quantum event, we consider two quantum observables Fq{^,x) and 
Gq{S,,x) and need to give meaning to the conditional expectation 

(5.1) E[F\G = g] 

We proceed in such a way that if the construction is apphed to ordinary probabiUstic conditioning 

and both the quantum density function and G arc diagonal, then tlic^ rc^sult is the same. Namely, 
we interpret the equation G = gwe start by an unconditioned density matrix pg and we construct 
the "closest" density matrix p| such that 

(5.2) TripSqG'^) = g" 

for all A; > 0. The point is to define what we mean by "closest". If for simplicity g is a discrete 
eigenvalue of G and Pg is the corresponding eigenspace, then a good definition is to set 



(5.3) p; 



g _ PgPqPg 

« rriPgPqPq)- 



This definition uniquely specifies Out of the quantum density matrix we can reconstruct 
many conditional classical densities. The lack of uniqueness of the construction is however 
immaterial since we are interested in taking expectations only of classical observables which are 
gauge invariant. 

All the standard quantum mechanics now follows. In particular, it is clear that two quan- 
tum observables F{$^,x) and G(^, x) are simultaneously measurable only if they commute when 
interpreted as matrices. 
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6. Conclusions 

We have defined a representation for non-relativistic quantum mechanics which is entirely 
equivalent to the standard theory but which is expressed in terms of a classical diffusion. 
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